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ABSTRACT 

The  analysis  of  B aimer-dominated  optical  spectra  from  nonradiative  (adiabatic)  SNRs  has  shown  that  the  ratio 
of  the  electron  to  proton  temperature  at  the  blast  wave  is  close  to  unity  at  vs  S  400  km  s  but  declines  sharply 
down  to  the  minimum  value  of  mjmp  dictated  by  the  jump  conditions  at  shock  speeds  exceeding  2000  km  s  '. 
We  propose  a  physical  model  for  the  heating  of  electrons  and  ions  in  non-cosmic-ray-dominated,  strong  shocks 
(vs  >  400  km  s  ')  wherein  the  electrons  are  heated  by  lower  hybrid  waves  immediately  ahead  of  the  shock  front. 
These  waves  arise  naturally  from  the  cosmic  ray  pressure  gradient  upstream  from  the  shock.  Our  model  predicts  a 
nearly  constant  level  of  electron  heating  over  a  wide  range  of  shock  speeds,  producing  a  relationship  ( TeITp)0  oc 
v;2  (oc M~2)  that  is  fully  consistent  with  the  observations. 

Subject  headings:  cosmic  rays  —  ISM:  kinematics  and  dynamics  —  plasmas  —  shock  waves  — 
supernova  remnants 


1.  INTRODUCTION 

The  discovery  of  collisionless  shock  waves  in  the  solar  wind 
in  the  1960s  ushered  in  a  new  era  in  the  physics  of  space 
plasmas  (see  Tidman  &  Krall  1971,  Sagdeev  1979,  and  Kennel 
1985  for  reviews  and  references).  Due  to  the  low  density 
(n  £  1  cm-3)  of  the  interplanetary  medium,  the  jump  in  hy- 
drodynamical  quantities  is  produced  not  by  Coulomb  collisions 
but  by  collective  plasma  processes  such  as  electromagnetic 
waves  and  turbulence.  However,  despite  the  availability  of  ex¬ 
tensive  in  situ  observations  of  solar  wind  shocks  and  the  ex¬ 
penditure  of  considerable  effort  in  theoretical  modeling  of  these 
structures,  a  detailed  understanding  of  the  processes  at  the 
shock  transition  responsible  for  partitioning  the  shock  energy 
between  different  charged  particle  species  has  been  slow  to 
emerge.  The  problem  is  far  more  acute  for  interstellar  shocks, 
where  in  situ  measurements  of  the  shock  structure  are  un¬ 
available  and  the  very  high  Mach  numbers  (~30-200)  make 
numerical  simulations  of  these  structures  extremely  difficult  or 
impossible. 

The  optical  emission  generated  by  nonradiative  supernova 
remnants  (SNRs)  in  partially  neutral  gas  provides  a  valuable 
diagnostic  tool  for  probing  the  heating  processes  in  collisionless 
shocks.  The  optical  spectra  of  these  SNRs  (which  lose  a  neg¬ 
ligible  fraction  of  their  energy  to  radiation)  are  dominated  by 
Balmer  line  emission,  produced  by  collisional  excitation  when 
neutral  hydrogen  is  overrun  by  the  blast  wave  (Chevalier  & 
Raymond  1978;  Bychkov  &  Lebedev  1979).  Each  emission 
line  consists  of  two  components:  (1)  a  narrow  velocity  com¬ 
ponent  produced  when  cold,  ambient  H  i  overrun  by  the  shock 
is  excited  by  electron  and  proton  collisions  and  (2)  a  broad 
velocity  component  produced  when  fast  neutrals  created  by 
postshock  charge  exchange  are  collisionally  excited  (Chevalier 
et  al.  1980).  An  example  of  a  Balmer-dominated  shock  spec¬ 
trum  from  the  Galactic  SNR  RCW  86  is  shown  in  Figure  1. 
The  optical  emission  arises  in  a  very  thin  (SlO16  cm)  ionization 
zone,  thin  enough  so  that  the  protons  transformed  into  hot 
neutrals  have  had  little  time  to  equilibrate  with  electrons  and 
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other  ions.  This  makes  the  measured  width  of  the  broad  Balmer 
line  directly  proportional  to  the  proton  temperature  set  by  col¬ 
lisionless  heating  at  the  shock  front.  The  broad-to-narrow  flux 
ratio,  /„//„,  on  the  other  hand,  is  sensitive  to  both  the  degree 
of  electron-proton  temperature  equilibration  at  the  shock  front 
[i.e.,  ( TJTp)0 ]  and  the  shock  velocity,  vs  (Chevalier  et  al.  1980; 
Smith  et  al.  1991).  The  ratio  also  depends  (although  less  sen¬ 
sitively)  on  the  preshock  neutral  fraction.  The  broad  component 
width  and  IB/IN  of  an  observed  Balmer-dominated  shock  can 
be  modeled  with  numerical  shock  codes  to  simultaneously  es¬ 
timate  vs  and  ( TJTp)0  (Smith  et  al.  1991;  Ghavamian  1999; 
Ghavamian  et  al.  2001,  2002). 

In  this  Letter,  we  draw  together  observed  values  of 
(' TJTp)0  and  vs  measured  in  five  Balmer-dominated  SNRs,  in¬ 
cluding  several  previously  unpublished  measurements  from  the 
SNR  RCW  86.  We  then  propose  a  physical  model  of  electron 
heating  by  lower  hybrid  waves  in  a  cosmic  ray  precursor  that 
obeys  the  observed  relationship  between  (TJTp)0  and  vs.  We 
conclude  by  exploring  some  consequences  of  this  interpreta¬ 
tion,  the  connection  to  other  observables,  and  applicability  to 
other  collisionless  shock  situations. 

2.  RELATIONSHIP  BETWEEN  EQUILIBRATION  AND  SHOCK  SPEED 

In  Figure  2,  we  plot  11  measurements  of  (TJTp)0  and  vs, 
obtained  from  long-slit  spectra  of  four  Galactic  remnants:  the 
Cygnus  Loop  (one  position  from  the  northeast),  RCW  86  (po¬ 
sitions  from  the  southwestern,  northern,  and  eastern  limbs),  Ty¬ 
cho’s  SNR  (“knot  g“  from  the  eastern  rim),  SN  1006  (one  po¬ 
sition  from  the  northwestern  rim),  and  one  remnant  in  the  Large 
Magellanic  Cloud,  DEM  L71  (five  positions  from  the  full  cir¬ 
cumference  of  the  shell).  In  the  first  four  remnants,  the  shock 
parameters  have  been  estimated  via  long-slit  spectroscopy.  First, 
we  used  the  broad  component  Ha  width  to  constrain  the  range 
of  shock  speeds  for  each  Balmer  filament  in  the  limits  of  minimal 
[{TJTp)0  =  me/mp ]  and  full  [{TJTp)0  =  1]  equilibration.  Next, 
we  fine-tuned  (TJTp)0  and  vs  by  using  shock  models  to  match 
the  observed  broad-to-narrow  flux  ratio  (Ghavamian  et  al.  2001, 
2002). 

In  the  case  of  RCW  86,  we  have  added  two  more  data  points, 
indicated  by  dashes  in  Figure  2,  derived  from  previously  un¬ 
published  optical  spectra  of  Ghavamian  (1999).  These  were  ob¬ 
tained  from  the  eastern  rim  [fFWHM^Ha)  =  640  ±  35  km  s-1; 
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Fig.  1. — Example  of  the  optical  spectrum  of  a  B aimer-dominated  shock, 
showing  the  broad  and  narrow  Ho:  lines  characteristic  of  nonradiative  shocks 
in  partially  neutral  gas.  This  spectrum,  originally  presented  by  Sollerman  et 
al.  (2003),  was  obtained  from  the  southwestern  rim  of  the  Galactic  SNR  RCW 
86,  with  high  enough  spectral  resolution  (~10  km  s_1)  to  resolve  the  broad 
(~500  km  s-1  FWHM)  and  narrow  (~30  km  s~‘  FWHM)  Ha  lines.  The  night- 
sky  OH  lines  (indicated  by  the  circled  plus  signs)  have  been  left  in  to  dem¬ 
onstrate  their  relatively  narrower  widths  compared  to  the  Ha  lines.  The  broad 
Ha  width  and  ratio  of  the  broad  to  narrow  Ha  flux  for  these  types  of  shocks 
were  used  to  produce  the  relationship  shown  in  Fig.  2. 

IB/IN  —  1.0  ±  0.2]  and  the  northern  rim  ['Gwiim^Hw)  = 
325  ±  10  km  s_1;  IB/IN  =  1.06  ±  0.1]  of  RCW  86.  In  the  fifth 
SNR  of  our  sample,  DEM  L71  (Ghavamian  et  al.  2003),  we 
used  the  broad  component  Ha  width  only  to  constrain  the  range 
of  shock  speeds  and  the  proton  temperature  Tp,  due  to  the  anom¬ 
alously  low  IB/IN  values  in  this  SNR  (see  below).  We  then  com¬ 
bined  this  information  with  Te  measured  from  Chandra  X-ray 
spectra  of  the  blast  wave  to  obtain  ( TJTp)0  (Rakowski  et  al. 
2003;  Rakowski  2005). 

Although  the  equilibrations  were  obtained  from  different  SNRs 
in  a  range  of  environments  and  distances,  they  show  a  clear  trend 
of  decreasing  ( TJTp)0  with  vs.  The  plot  in  Figure  2  suggests  that 
a  transitional  shock  speed  exists  around  400  km  s~\  below  which 
collisionless  processes  promptly  equilibrate  the  electrons  and  pro¬ 
tons  at  the  shock  front.  Above  that  speed,  (TJTp)0  rapidly  declines, 
eventually  reaching  values  consistent  with  mass-proportional  heat¬ 
ing.  This  decline  appears  to  follow  ( TJTp)0  oc  vs2 .  We  show  that 
this  behavior  is  predicted  by  a  physical  model  of  electron  heating, 
where  the  electrons  are  heated  in  the  cosmic  ray  precursor  to  a 
level  that  is  nearly  independent  of  shock  speed  before  acquiring 
the  mass-proportional  increment  at  the  shock  front.  The  protons, 
on  the  other  hand,  receive  only  mass-proportional  heating,  re¬ 
sulting  in  the  inverse  squared  dependence  of  equilibration  on  shock 
speed.  Note  that  while  our  proposed  model  relies  on  the  accel¬ 
eration  of  cosmic  rays  and  the  existence  of  a  cosmic  ray  precursor 
to  heat  the  electrons,  it  does  not  require  the  energetics  of  the  shock 
to  be  dominated  by  the  cosmic  rays. 

In  physical  models,  the  strength  of  the  (assumed  quasi- 
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Fig.  2. — Electron  to  proton  temperature  ratio  at  the  shock  front  as  a  function 
of  shock  velocity  for  five  B aimer-dominated  SNRs.  Magnetosonic  Mach  numbers 
(Ms)  appropriate  for  typical  ISM  conditions  are  indicated  along  the  top  axis.  The 
data  shown  here  were  measured  from  Balmer-dominated  shocks  in  the  Cygnus 
Loop,  RCW  86,  Tycho’s  SNR  (Ghavamian  et  al.  2001),  SN  1006  (Ghavamian 
et  al.  2002),  and  DEM  L71  (Rakowski  et  al.  2003;  Rakowski  2005).  The  dashed 
error  bars  for  RCW  86  mark  previously  unpublished  results.  Below  400  km  s-1 
(Ms  «  30),  the  data  are  consistent  with  (Te/Tp)0  =  1.  The  prediction  of  the  pro¬ 
posed  lower  hybrid  wave-heating  mechanism  in  the  cosmic  ray  precursor, 
(Tj/T,,) „  oc  V“2  (ocMj-2),  is  shown  for  vs  >  400  km  s_l. 

perpendicular)  shock  is  characterized  not  by  the  shock  speed  but 
rather  by  the  magnetosonic  Mach  number  Ms  (=vs/vMS,  where 
vMS  =  (cB  +  vA)ln  is  the  magnetosonic  speed,  cs  is  the  sound 
speed  (  =  [(5/3)(R/p)]1/2)  and  vA  [  =B/(47rp,)1/2]  is  the  Alfven 
speed  of  the  preshock  gas).  The  preshock  temperature,  ion  den¬ 
sity,  and  magnetic  field  strength  are  not  strongly  constrained  in 
the  observed  Balmer-dominated  shocks.  In  particular,  Ms  is  most 
sensitive  to  the  choice  of  preshock  magnetic  field  due  to  the 
dependence  of  the  Alfven  speed  onB2.  However,  if  we  assume 
standard  values  for  the  warm  neutral  interstellar  medium 
(ISM) — preshock  temperature  of  10,000  K,  density  of  1  cm-3, 
magnetic  field  strength  of  3  /xG,  and  50%  preshock  ionization — 
the  magnetosonic  speed  is  then  approximately  13  km  s_1.  The 
corresponding  values  of  Ms  are  marked  at  the  top  of  Figure  2. 

3.  LOWER  HYBRID  WAVE  HEATING  MODEL 

The  constant  electron  heating  with  shock  velocity  [giving 
(TJT  ) „  oc  vy2~\  suggests  a  process  occurring  in  the  preshock 
medium  rather  than  the  shock  front  itself.  Waves  in  the  pre¬ 
shock  medium  can  be  excited  by  shock-reflected  ions,  which 
gyrate  around  the  field  lines  before  returning  to  the  shock  itself, 
or  by  a  cosmic  ray  precursor.  Cargill  &  Papadopoulos  (1988) 
performed  hybrid  simulations  of  the  first  possibility,  wherein 
shock-reflected  ions  generate  Langmuir  and  ion-acoustic  waves 
ahead  of  the  shock,  which  then  heat  electrons  as  they  are 
damped.  They  predicted  a  temperature  ratio  (TJTp)0  ~  0.2 
nearly  independent  of  shock  velocity,  a  result  that  clearly  dis¬ 
agrees  with  the  observational  data  in  Figure  2.  Laming  (2001a, 
2001b)  modeled  the  excitation  of  lower  hybrid  waves  by  shock- 
reflected  ions,  following  the  suggestion  of  McClements  et  al. 
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(1997)  that  lower  hybrid  waves  could  stream  away  from  the 
shock  with  a  group  velocity  equal  to  the  shock  velocity.  In  this 
way,  waves  stay  in  contact  with  the  shock  for  arbitrarily  long 
periods  of  time  and  hence  can  grow  to  large  intensities  despite 
low  intrinsic  growth  rates.  The  electron  heating  predicted  by 
such  a  model  (Vink  &  Laming  2003)  depends  on  the  product 
of  the  maximum  electron  energy  and  the  fraction  of  electrons 
that  are  resonant  with  the  waves  that  can  be  accelerated.  This 
gives  (Te/Tp)0  oc  v j 1  exp  (—  M2),  where  M  is  the  shock  Mach 
number. 

Here  we  explore  another  approach.  Cosmic  rays  have  been 
long  known  to  generate  waves  upstream  of  shocks;  the  gen¬ 
eration  of  Alfven  waves  is  an  intrinsic  part  of  cosmic  ray  ac¬ 
celeration  models.  Drury  &  Falle  (1986)  showed  that  the  neg¬ 
ative  gradient  of  cosmic  pressure  with  distance  ahead  of  the 
shock  can  generate  sound  waves  via  the  ''Drury  instability.” 
This  is  the  mechanism  that  smooths  out  the  hydrodynamical 
jump  in  cosmic  ray  modified  shocks.  Mikhailovskii  (1992) 
gives  a  corresponding  expression  for  the  excitation  of  mag¬ 
netoacoustic  waves.  We  argue  that  the  high-frequency  exten¬ 
sion  of  such  processes  by  the  cosmic  ray  pressure  gradient 
would  be  the  excitation  of  lower  hybrid  waves.  The  frequencies 
of  these  waves  would  lie  between  the  gyrofrequencies  of  the 
electrons  and  the  protons,  with  electron  heating  occurring  as 
the  waves  damp  along  magnetic  field  lines.  Without  calculating 
the  growth  rate  of  the  lower  hybrid  waves  explicitly,  we  can 
estimate  the  magnitude  of  the  electron  heating.  The  parallel 
diffusion  coefficient  for  electrons  in  lower  hybrid  wave  tur¬ 
bulence  is  (derived  from  eqs.  [10.83]  and  [10.93]  of  Melrose 
1986;  see  also  Begelman  &  Chiueh  1988) 


where  co  =  (O^p)1^  is  the  lower  hybrid  wave  frequency,  the 
geometric  mean  of  the  electron  and  proton  cyclotron  frequen¬ 
cies  (Q„  p  =  eB/me  pc),  and  fey  and  kL  are  wavevectors  parallel 
and  perpendicular  to  the  magnetic  field.  Taking  8E  = 
B  (fip/co)1/3  (co/4 fcjc)  (Karney  1978),  =  melmp,  and  the 

time  spent  by  an  electron  inside  the  cosmic  ray  precursor  t  ~ 
k!v%  where  k  is  the  cosmic  ray  diffusion  coefficient,  we  find 


would  be  t  ~  dlvs  ~  1/Hp,  where  the  precursor  depth  d  is  ap¬ 
proximately  the  gyroradius  of  a  cosmic  ray  proton.  Substituting 
this  relation  into  equation  (2)  gives  ( TJTp)0 ,  independent  of 
vs,  similar  to  the  behavior  predicted  by  Cargill  &  Papadopoulos 
(1988),  although  these  authors  consider  different  waves.  In  this 
case,  the  shallowness  of  the  reflected  ion  precursor  prevents 
the  collisional  redistribution  of  energy  among  the  electrons. 

If  a  constant  kinetic  energy  is  imparted  to  electrons  by  lower 
hybrid  waves,  then  for  shock  speeds  of  400  km  s  1  and  above 
(cf.  Fig.  2)  the  observations  require  mev2J2  ~  4  x  1CT10  ergs 
(~0.3  keV)  immediately  behind  the  shock  front.  In  that  case, 
we  infer  UeK  -  7  x  1021  cm2  s~2,  and  k  ~  4  x  1014/B  cm2  s_1. 
For  typical  magnetic  fields  B  ~  3  /xG,  k  is  several  orders  of 
magnitude  below  that  inferred  for  the  undisturbed  ISM  but 
entirely  consistent  with  values  estimated  for  the  solar  wind  or 
interplanetary  medium.  We  emphasize  that  this  estimate  is  a 
lower  limit  on  k.  This  estimate  neglects  the  electron-electron 
collisional  equilibration  necessary  to  bring  more  electrons  into 
resonance  with  the  turbulence.  Furthermore,  lower  hybrid 
waves  can  only  be  excited  in  the  cosmic  ray  precursor  if  the 
cosmic  ray  pressure  gradient  is  sufficiently  strong.  The  required 
gradient  scale  lengths  are  1/L  >  8c5/3k  for  sound  waves  (Drury 
&  Falle  1986)  and  1 IL  >  vA/crg  ~  vJk  for  magnetoacoustic 
waves  (Mikhailovskii  1992),  where  rg  is  the  gyroradius  of  a 
cosmic  ray  proton.  An  upper  limit  on  k  will  come  from  the 
requirement  that  the  neutral  H  survive  against  electron  impact 
ionization  in  the  precursor  to  reach  the  shock  front.  Numeri¬ 
cally,  xS  1024  (i^/lOOO  km  s^1)2 /ne  cm2  s~\  Further  con¬ 
straints  on  the  cosmic  ray  diffusion  coefficient  and  pressure 
will  be  obtained  from  calculation  of  the  lower  hybrid  wave 
growth  rate  required  to  sustain  the  electron  heating.  This  growth 
rate  must  be  significantly  larger  than  the  charge  exchange  fre¬ 
quency  of  the  partially  neutral  gas.  Under  typical  ISM  con¬ 
ditions,  (f^flp)172  ~  1  Hz  and  nHo(acxv)  ~  1CT8  Hz,  so  the  con¬ 
dition  for  electron  heating  is  nearly  always  satisfied. 

The  electron  heating  by  lower  hybrid  waves  is  inherently 
anisotropic.  If  some  of  this  anisotropy  survives  the  electron- 
electron  collision  equilibration,  a  polarization  signal  may  be 
present  in  the  narrow  component  of  Ha  (see,  e.g..  Laming 
1990),  which  might  provide  more  insight  into  the  electron  heat¬ 
ing  process. 
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If  the  lower  hybrid  wave  group  velocity  dw/c )kL  =  vs,  then  the 
phase  velocity  oi/k±  =  2vs  when  kpk]_  =  melmp  (Laming 
2001a),  and  the  electron  heating  depends  on  a  few  constants 
times  the  product  QeK.  Thus,  for  a  Bohm-like  cosmic  ray  dif¬ 
fusion  coefficient  k  oc  MB  and  independent  of  vs,  we  have  an 
electron  heating  process  independent  of  both  the  shock  velocity 
and  the  upstream  magnetic  field,  and  ( TJTp)0  oc  vs  1  as  required. 
This  estimate  neglects  the  fact  that  only  electrons  with  velocity 
greater  than  ve  ~  u/k^  may  be  heated  by  the  waves.  However, 
for  the  likely  depth  of  a  cosmic  ray  precursor  (estimated  below), 
the  electrons  will  have  sufficient  time  to  at  least  partially  col- 
lisionally  equilibrate  among  themselves  before  crossing  the 
shock,  allowing  a  much  larger  fraction  of  preshock  electrons 
to  interact  with  the  waves.  The  time  spent  by  an  electron  inside 
a  reflected  ion  precursor,  as  opposed  to  a  cosmic  ray  precursor. 


4.  THE  WIDTH  OF  THE  Ha  NARROW  COMPONENT 

The  same  cosmic  ray  precursor  should  also  generate  lower 
frequency  waves,  which  for  a  quasi-perpendicular  shock  will 
include  magnetoacoustic  waves.  Below  the  ion-neutral  charge 
exchange  frequency,  these  waves  are  not  effectively  damped 
(Drury  et  al.  1996)  and  may  reveal  themselves  via  broadening 
of  the  narrow  Ha  component  (Smith  et  al.  1994).  Those  authors 
suggested  that  waves  in  the  cosmic  ray  precursor  actually  heat 
the  preshock  gas  to  30,000-40,000  K  and  that  the  observed  line 
width  is  thermal  in  nature.  An  upper  limit  to  the  cosmic  ray 
diffusion  coefficient  of  ~1024/ne  cm2  s  then  results  from  the 
constraint  that  sufficient  neutral  H  survive  passage  through  the 
precursor  to  encounter  the  shock.  However,  no  precise  heating 
mechanism  was  specified,  and  we  now  suggest  that  the  large 
width  of  the  narrow  component  Ha  line  in  the  observed  SNRs 
is  not  of  thermal  origin  but  rather  due  to  the  motion  of  protons 
in  the  lowest  frequency  waves  of  the  magnetoacoustic  spectrum. 
These  waves  lie  below  the  charge  exchange  frequency,  allowing 
a  coherent  oscillation  of  preshock  neutrals  and  preshock  protons. 
In  this  case,  8v  —  vA8B/B  for  Alfven  waves  (i>A  should  be  replaced 
by  vMS  for  magnetoacoustic  waves).  The  Alfven  speed  upstream 
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from  the  shock  is  in  the  range  1-10  km  s  so  an  observed 
broadening  8v/vA  >  1  also  implies  8B/B  >  1 .  Such  a  magnetic 
field  amplification  has  already  been  proposed  in  the  context  of 
cosmic  ray  acceleration  in  shocks  (Lucek  &  Bell  2000;  Bell  & 
Lucek  2001).  A  successful  model  of  precursor  ions  will  need  to 
reproduce  the  existing  measurements  of  narrow-component  Ha 
line  widths  (Smith  et  al.  1994;  Hester  et  al.  1994;  Sollerman  et 
al.  2003),  which  indicate  that  they  remain  relatively  constant 
(30-50  km  s  ')  over  a  wide  range  in  shock  speed  (180-3000 
km  s  ').  We  defer  this  topic  to  future  work. 

5.  DISCUSSION  AND  SUMMARY 

The  plot  of  (TJTp)0  in  Figure  2  is  qualitatively  similar  to 
results  of  a  survey  of  electron  heating  at  solar  wind  shocks 
(Schwartz  et  al.  1988).  The  main  difference  is  that  Te  ~  7}  is 
obtained  for  Alfvenic  Mach  numbers  <5  in  the  solar  wind, 
whereas  in  Figure  2  this  occurs  for  magnetosonic  Mach  numbers 
around  20-30,  for  assumed  preshock  magnetic  fields  of  3  /xG. 
A  preshock  magnetic  field  amplification  in  SNRs,  of  a  similar 
amount  required  to  reproduce  the  width  of  the  Ha  narrow  com¬ 
ponent  (8B/B  ~  5-10),  would  shift  the  ( TJTp)0  ratios  in  Fig¬ 
ure  2  onto  Mach  numbers  similar  to  those  in  the  solar  wind. 
This  may  support  our  interpretation.  Collisionless  shocks  play  a 
dominant  role  in  heating  the  intergalactic  medium  (IGM)  during 
large-scale  structure  formation  (Ryu  et  al.  2003;  Yoshida  et  al. 
2005;  Kang  et  al.  2005).  If  the  preshock  magnetic  field  is  am¬ 
plified  as  described  above,  then  the  relationship  ( TJTp)0  oc 
Ms  2  may  also  be  applicable  to  modeling  the  ionization  and  tem¬ 
perature  structure  of  the  IGM. 

The  maximum  energy  attainable  by  cosmic  rays  in  B  aimer- 


dominated  shocks  is  limited  by  ion-neutral  damping  to  ~0.1- 
1  TeV  (Drury  et  al.  1996).  It  is  possible  that  the  shocks  with 
stronger  cosmic  ray  acceleration  generate  sufficient  electron 
heating  in  their  precursors  that  no  neutrals  survive  to  cross  the 
shock.  The  anomalously  low  IB/IN  ratio  observed  in  DEM  L71 
(Ghavamian  et  al.  2003)  may  indicate  that  this  is  an  interme¬ 
diate  case.  Electron  heating  in  a  shock  precursor  will  produce 
added  narrow-component  Ha.  However,  since  the  broad  com¬ 
ponent  can  only  arise  by  charge  exchange  with  shocked  pro¬ 
tons,  the  neutral  H  must  penetrate  either  the  shock  or  the  re¬ 
flected  ion  precursor.  We  emphasize  that  of  the  points  plotted 
in  Figure  2  from  Ha  /„//v  measurements,  the  electron  tem¬ 
peratures  from  RCW  86  (Vink  et  al.  2006),  the  Cygnus  Loop 
(Levenson  et  al.  2002),  SN  1006  (Laming  et  al.  1996;  Vink  et 
al.  2003),  and  Tycho  (Hwang  et  al.  2002;  Warren  et  al.  2005) 
have  been  independently  corroborated  by  measurements  in 
other  wave  bands. 

Several  lingering  questions  remain  from  Figure  2:  What 
mechanism  causes  the  prompt  electron-proton  equilibration  be¬ 
low  400  km  s-1,  and  why  does  the  decline  in  ( TJT  )0  begin 
at  that  shock  speed?  Even  more  importantly,  does  the  inverse 
square  relationship  between  equilibration  and  shock  speed/ 
Mach  number  also  hold  for  collisionless  shocks  in  fully  ionized 
gas?  These  issues  will  be  addressed  in  future  work. 
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